Evaporative cooling in a radio— frequency trap 
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A theoretical investigation for implementing a scheme of forced evaporative cooling in radio- 
frequency (rf) adiabatic potentials is presented. Supposing the atoms to be trapped by a rf field 
at frequency wi, the cooling procedure is facilitated using a second rf source at frequency L02 ■ This 
second rf field produces a controlled coupling between the spin states dressed by U\ . The evaporation 
is then possible in a pulsed or continuous mode. In the pulsed case, atoms with a given energy are 
transferred into untrapped dressed states by abruptly switching off the interaction. In the continuous 
case, it is possible for energetic atoms to adiabatically follow the doubly-dressed states and escape 
out of the trap. Our results also show that when uji and UJ2 are separated by at least the Rabi 
frequency associated with ui, additional evaporation zones appear which can make this process 
more efficient. 

PACS numbers: 32.80.-t, 39.25.+k, 32.80.Pj 



I. INTRODUCTION 

In recent years, the investigation of quantum gases 
in low-dimensional trapping geometries has significantly 
attracted the attention of the physics research commu- 
nity. This growing interest is motivated, partially, by the 
current possibilities that the extremely rapid progress 
in integrated atom optics has opened for the manipu- 
lation of Bose-Einstein condensed (BEC) atoms. This 
development allows the study of crucial problems asso- 
ciated with the strong modifications that the fundamen- 
tal properties of these quantum systems experience due 
to the reduced dimensionality. For instance, a ID Bose 
gas in the Tonks-Girardeau regime mimics a system of 
non-interacting spinless fermions 0,0,0]; in 2D, the su- 
perfluidity emerges due to the vortex binding-unbinding 
Berezinskii-Kosterlitz-Thouless phase transition 0, 
recently observed [(|. 

For the study of the BEC low-dimensional physics, 
trapping configurations of different nature and topology 
have been proposed and used. For example, the 3D to 
ID crossover was explored by Gorlitz et al. |jj in an elon- 
gated Ioffe-Pritchard type direct-current (dc) magnetic 
trap, the phase defects of a BEC were investigated in 
a quasi-2D trap based on a ID optical lattice [8( and, 
in atom chip experiments, dc current-carrying wires are 
usually employed to confine atoms in highly anisotropic 
traps Although these trapping configurations have 
demonstrated their relevance for studying quantum gases 
in low dimensions, adiabatic potentials jlO(, resulting 
from a combination of dc and radio-frequency (rf) mag- 
netic fields, are also becoming a very attractive and 
promising tool El El El El 111 E3 

The rf traps share the versatility and flexibility of the 
above mentioned trapping schemes and, moreover, they 



'Electronic address: leonardo@galilce.univ-parisl3.fr 



are relatively easy to implement and control. In the first 
implementation of these traps , ultra-cold atoms were 
confined in a 2D geometry. A rf adiabatic potential has 
also been used as a beam splitter, allowing the demon- 
stration of matter-wave interference on an atom chip |l2l | . 
Ring-shaped traps, and other more complex trapping ge- 
ometries using adiabatic potentials have also been con- 
sidered [Hlil El [3. 

Given the topology of the rf trapping potential, and be- 
cause of technical limitations in some cases, the loading of 
the trap with Bose-Einstein condensed atoms, preserving 
the quantum degeneracy, can be a challenging task. In 
this situation, it is of relevance to consider the possibil- 
ity of evaporative cooling of atoms directly in these low- 
dimensional rf traps. This is the subject that will be ad- 
dressed in this paper, taking into account the interaction 
of the atoms with two radio-frequency fields. When deal- 
ing with more than one rf frequency, an analytical solu- 
tion for the atomic spin dynamics can be found by treat- 
ing the individual successive interactions of the rf fields 
with the atoms 17] or by considering the two fields simul- 
taneously, provided one of the fields is rather weak jl8| . 
We will study how a weak second radio-frequency source 
can be used to perform an evaporation. 

We will see that the forced evaporation of rf-trapped 
atoms can be accomplished in two ways. Firstly, the 
spin evolution induced by this second rf source can be 
quenched by switching off the field, i. e. by using a pulsed 
rf source. Secondly, we can allow an adiabatic follow- 
ing of doubly-dressed states which requires the second rf 
source to be continuous rather than pulsed. This last 
scheme is similar to the standard evaporative cooling 
method used in static magnetic traps El LSI \2M ■ 

This paper is organized as follows: In section [H] we 
will discuss the geometry of the system and the singly- 
dressed states of the rf trap. In Sec. IIIII the evolution 
of the system is determined in three different ways: nu- 
merically, using a first order Magnus series approxima- 
tion, and by using a second rotating wave approximation 
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which leads to a double-dressing of the atoms by two rf 
fields. Sec. I1VI is devoted to the application of the re- 
sults of Sec. lIIII to the study of evaporative cooling in the 
rf trap. Finally, we give a summary and conclusion in 
Sec. E| 

II. ADIABATIC POTENTIAL CONFINEMENT 

The underlying idea of the confinement of ultra-cold 
atoms using rf adiabatic potentials is presented in detail 
in Ref. 01 > ^ er - E3 being the first report on the experi- 
mental investigation of such a trapping scheme. For this 
reason, instead of discussing deeply how this trapping ac- 
tually takes place, we will rather make use of the already 
known results that are relevant in order to consider the 
problem of evaporative cooling in these traps. 

The treatment presented in this paper is valid for any 
value of the spin F. However, for concreteness the nu- 
merical results will be given for 87 Rb in the Zeeman state 
rn = +2 of the 551/2, F = 2 hyperfine ground state level. 
We will suppose that the atoms are confined in a QUIC 
magnetic trap (2^, produced by a dc magnetic field 
Bd c (r). The atomic clouds trapped in this configuration 
are anisotropic (cigar-shaped along x) and we will take 
the offset magnetic field produced by the Ioffe coil 22] to 
be oriented along the x direction. In the following, the 
axes in the lab frame will be labelled by lower case letters 
xyz. The z axis is in the vertical direction and y is the 
horizontal direction perpendicular to the cigar axis. The 
axes in the local frame attached to the dc magnetic field 
will be labelled by capital letters XYZ. Moreover, we 
will assume that the direction of the dc magnetic field 
defines the local Z quantization axis. The Larmor fre- 
quency of the atomic spin precession in this dc field will 
be denoted loq(y) = gF^BBdc^) /h. Here, gF and /is are 
the Lande factor and the Bohr magneton, respectively. 

Now, we apply to this confined atomic system two rf 
fields (produced by antennae), both of them polarized 
along X and of angular frequencies lo\ and uj2- Having in 
mind that the second rf field will be rather weaker than 
the first one, we transform into the frame rotating at oj\ 
and perform the rotating wave approximation (RWA). 
The Hamiltonian that describes the spin dynamics can 
thus be written as (see the derivation in appendix |A"|) 

H{r, t) = H A (r) + Q 2 [F x cos(Ai) + F Y sin(At)] , (1) 

where Fx(Fy) is the atomic angular momentum in the 
X(Y) direction, A = oj 2 - 0Ji, H A (r) = tt(r)F g is 
the adiabatic Hamiltonian associated with the rf con- 
finement, and ^(r) = ^/S(r) 2 + Vl\ defines the energy 
separation between the adiabatic levels. In the ab- 
sence of O2 , the flip angle 9 and the detuning 5(r) 
are given by tan(#) = — Oi/ <5 (r) , with 9 6 [0, 7r], and 
8(r) = u>i — u>o(t), respectively. We have labelled the 
Rabi frequencies of the rf fields Sli and SI2, and we have 
considered that the component Fz of the atomic angular 
momentum is aligned with the local Z component of the 



dc magnetic field vector. Strictly, these Rabi frequen- 
cies are not spatially homogeneous, however they can be 
treated as such over the spatial extension of the atomic 
cloud [13 • 

Graphically, the spin evolution given by Eq. Q), in the 
case f^2 = 0, is represented in Fig. It can be seen in 
this figure that the tilted angular momentum F$ results 
from a rotation of Fz around Fy and is given by 

F e = cos(9)F z + sm(e)F x = K Y {9)F Z TZ Y (9) , (2) 

where the rotation matrix TZy{9) = exp(—i9Fy/h) can 
be expressed in the basis {—2, +2} of the bare states. 




FIG. 1: (color online). At a given location r, S and Qi define 
the angle 9 by which Fz is rotated. The spin (black arrow) 
therefore precesses around an axis given by Fe (green arrow) 
at a frequency fi. 



In Fig. Eta), the energies of the bare states are plotted 
as a function of the position z, where the energy varia- 
tion is due to the dc magnetic field B<j c (r). This spatial 
dependence has been calculated for a value x = x m in = 
6.9 mm corresponding to the position where the QUIC 
magnetic field is minimal in our experimental setup . 
Moreover, we have taken y — and a Rabi frequency 
f2i/27r = 400 kHz. The arrows shown in that figure, in 
blue for lj\ and in red for ix>2 (each frequency stands to 
the right of its respective arrows), indicate the locations 
where the corresponding rf fields resonantly couple the 
states in the laboratory frame. On the other hand, the 
spatial z dependence of the adiabatic states internal en- 
ergies is shown in Fig.EJb) where the states are labelled, 
from top to bottom: |+2^), | + \Qa), | — 1a), and 
I — 2a) ■ We can also see in this last figure the avoided 
level crossings at the positions where u)\ (taken equal to 
27rx3.19 MHz in this example and from now on) reso- 
nantly couples the bare states. 

In order to consider only the confinement (f^ = 0) 
using the adiabatic potentials shown in Fig. Hfb), let's 
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FIG. 2: (color online). Spatial z dependence of the energy of 
the uncoupled (a) and dressed (b) states for x = x m in and 
y = 0. At the avoided level crossings, the energy splitting 
between the levels is fii. Gravity has been taken into account 
in both cases and LJ2 > u>j . 



suppose that initially we have am — +2 spin polarized 
ultra-cold atomic sample. In this situation, the trapping 
potential corresponding to the bare state | + 2) can be 
adiabatically deformed into the one associated with the 
dressed state | + 2 a)- Such a transformation can be per- 
formed by changing the detu ning 8( r) from red to blue at 
constant Rabi frequency f2i llfilll| or, by increasing Qi 
at a constant red detuning [l2|. Here, by adiabatic de- 
formation we mean that the angular precession frequency 
f2(r) of the spin in Fig.^niust be much larger than the 
rate at which the angle 8 changes [\9\ <C f2(r)]. Using the 
loading schemes just mentioned, it is possible to obtain 
highly anisotropic rf traps with trapping frequencies, in 
the strongest confinement direction, ranging from several 
hundreds of Hz up to a few kHz. 

Having described the main properties of the adiabatic 
confinement, let's now address the following issues. As- 
suming intuitively the existence of the resonances repre- 
sented by the arrows in Fig.^b), we would like to know 
precisely where they are located and how strong they are. 
Another relevant point to be taken into account concerns 
the effect of these resonances at the rf trap centre when, 
numerically, u>x and u>2 are close to each other. Moreover, 
it will be interesting to find out the different parameter 
values for which the second rf field induces transitions 
between the adiabatic states, in the perturbative limit 
with f^2 <C leading to a possible implementation of 
evaporative cooling in rf traps. 



III. DYNAMICS OF THE SYSTEM 

In this section we will study the dynamics of the system 
using three different methods. In the first case, the exact 
numerical solution of the time-dependent Schrodinger 
equation (TDSE) will be found. Secondly, an approx- 
imated analytical treatment will be presented (Magnus 
approximation) in order to interprete the exact results 
derived numerically. Finally, we will present an ana- 
lytic solution based on a second rotating wave approx- 
imation which will be the basis of the analysis presented 



in SeclTVl 

Since the evolution of the atomic external and inter- 
nal degrees of freedom takes place on very different time 
scales, here we will decouple the two dynamics and con- 
sider only the time evolution of the internal degrees of 
freedom. 



A. Numerical solution 

The evolution of the state vector ^(t)) with the Hamil- 
tonian was solved numerically, in the interaction pic- 
ture, using a 4th order Runge-Kutta algorithm. In this 
case, the state vector can be very efficiently propagated 
in time and we do not expect to have convergency prob- 
lems if we choose an appropriate time-step |24| . The first 
question we would like to address here is: supposing that 
an atom is initially in the trapped dressed state | +2^}, 
what is the probability Pi a = (2 a ^(r, t))\ 2 of finding it 
in that same state as time goes by? We will also be inter- 
ested in how this probability changes for an atom located 
in different places in the trap. The preliminary answer 
to these questions is presented in Fig. El where the prob- 
ability we are interested in is plotted for three different 
values of A, the difference between the radio-frequencies 
lo\ and UJ2- 

In Fig. El and from now on, we use fix as a frequency 
unit, having in mind to keep it constant in a given exper- 
imental situation. The value of has been taken equal 
to 0.05f2i, small enough in order to be in the perturba- 
tive limit, as will be shown later. This value of O2 will 
be used in all the following results unless a different one 
is explicitly stated. In the left column of Fig. El the x 
and time dependence of P2A are calculated at the points 
y = and z m i n k, —0.19 mm, with this value of z corre- 
sponding to the location of the left avoided level crossing 
in Fig. |2fb) - To obtain the right column of figure 3, we 
use for x the value x m i n ~ 6.9 mm, very close to the 
location where the QUIC trap has its minimum in the x 
direction. 

In Fig. El the three values of A have been chosen to 
illustrate some key characteristics of the spin dynamics. 
As the energy separation between the adiabatic levels at 
the avoided level crossings is exactly fix, we observe res- 
onant behaviour at the rf trap bottom for 0J2 = u>i + fii, 
i.e. when A = fii [see Fig. Etc) and (d)]. Similarly, for 
A < ill we have a red-detuned interaction everywhere 
[Fig. Efa) and (b)]. In this case we observe weak modu- 
lations of P2A that are essentially determined by a beat- 
ing between the A and f2(r) frequency components. For 
A > Hi [Fig-Et e ) an d (f)] we have a blue-detuned inter- 
action around the minima x m i n and z m i n . Away from the 
centre a resonance occurs, as expected, at approximately 
the location in the trap where 0J2 resonantly couples the 
bare states [see Fig. Et a )]- This is the outer resonance 
labelled OR in Fig. Eff). However, there is another fea- 
ture we would like to stress. Namely, the presence of the 
inner resonance IR clearly seen in the z dependence of 
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FIG. 3: (color online). Probability for an atom to remain 
in the initial rf trapped state | + 2a). The values of A are 
respectively: 0.875Sli in (a) and (b); fii in (c) and (d); and 
1.25f2i in (e) and (f). The rf trap frequency u>i/2tt is set to 
3.19 MHz. The x-t (z-t) dependence of P2A is calculated at 
Zmin = —0.19 mm (xmin = 6.9 mm). In (f) the labels OR and 
IR indicate the locations of the outer and inner resonances, 
respectively. 



P2A in Fig- 130- Note that the avoided level crossing 
(rf trap centre) is at z m i n . Looking back to Fig. and 
having in mind that L02 > wi, the existence of this second 
resonance IR in Fig. Elf) at z > z m i n may be counter- 
intuitive and, as we will see, its relative strength is fully 
determined by the rotation angle 9. Note that because 
of the loose confinement in the x direction [Tlj . the dy- 
namics in the y-z plane does not change much from one 
location to another in the x axis. However, this dynam- 
ics is very sensitive to changes in z (or y) and therefore, 
the results for the x dependence of the probability Pia in 
Fig. Eta), (c), and (e) can be significantly different when 
another z location is considered. 



B. Interpreting the numerical results using a first 
order Magnus series approximation 

Searching for the understanding of the physical picture 
behind the numerical results presented in Fig. El let's con- 
sider the first order Magnus series approximation [2^ |2|j 
to the solution of the TDSE. This approximation is ba- 
sically the formal solution of the TDSE neglecting the 
two-time commutators of the Hamiltonian. This Hamil- 
tonian is given in the interaction picture by 

H (t) = exp{intF e /h)H(t) exp(-iQtF g /h) . (3) 

In Eq. J3J), we have dropped the r dependence in fi 
and Fg for the sake of notational simplification. Making 
use of Fig. [I] and introducing F±g = K Y {6)F x n Y {6) = 
cos(8)Fx — sin(8)Fz (see Fig. EJ, which is the angular 
momentum vector perpendicular to Fg and Fy, wc find 



H{t) = n 2 {sin(0) coa{M)F g + [sin(At) sin (fit) + cos(0) cos(Ai) cos(fit)]Fj_ e + 

[ sin(At) cos(fit) - cos(0) cos( At) sin (fit)] F Y | . (4) 
I 



In this case, the time evolution of a given initial dressed 
spin state |^(0)) is represented by the rotation 

|*(f))=exp(-iE(t).F/ft)|tt(0)> , (5) 

where the scalar product H(t) • F stands for Eg(t)Fg + 
E ± g(t)F ±e + E Y (t)F Y , where E g (t), E±g(t) and E Y (t) are 
a measure of the projections of the spin |\&(t)) precession 
axis onto the axes Fg, F±$, and F Y , respectively (for 



simplicity we will just call them projections) . Taking into 
account the equation lf4"|. the definition of the exponential 
argument in Eq. |JSJ is 



3(t)-F= / dt'H'(t') , (6) 
Jo 

leading to the following expressions for the projections 



Se(t) = sin(0) sin(At) , 



(7) 
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3±e 0) =ft 2 {cos 2 (0/2) 



sin[(A - Q)t] 

a - n 



-sin 2 ((9/2) 



sin[(A 



(8) 



E Y (t) = 2n 2 1 cos 2 (6»/2) 



sin 2 [(A - Q)t/2] 

a- n 



sin 2 (6>/2) 



sin 2 [(A + fi)t/2] 
A + ft 



(9) 



r 



By inspecting Eqs. (JZJh©, we can see the time- 
dependent terms resulting from a beating between fre- 
quency components at A and f2(r). These beats are seen 
as the modulation (interference-like patterns) of P%a ob- 
served in Fig. |21 In particular, we sec in these equa- 
tions that there will be some interesting behaviour when 
A = ±fi. In either case the condition is realized by two 
values of 8: 8q and tt — 8q with 9q = arcsin(Sli/A). If 
A ^ ±fi, the coefficients Bj(i) are oscillatory with finite 
amplitudes. However, when A = Q at 9 = 8q, for in- 
stance, Sj_e(t) shows a linear tendency in time of the form 
f^2 cos 2 (6 l o/2)t while the other two components are neg- 
ligible. This suggests that the spin will essentially rotate 
at a frequency fl 2 cos 2 (6*o/2) around the axis F±g. Start- 
ing from an eigenstate of Fg, as in Sec. IIII Al the spin 
will be completely flipped after a half period. This rota- 
tion corresponds to the outer resonance OR in Fig- Elf) - 
At the location 6 = n — 9q, the same resonant behaviour 
occurs with a rotation frequency sin 2 (6V)/2). This cor- 
responds to the inner resonance IR in Fig. EJ^f ) . In the 
case when 0J2 is smaller than uji, we have the resonant 
condition A = — fi and we have the same behaviour ex- 
cept that the character of the inner and outer resonances 
is now reversed. 

Away from the resonant conditions just described the 
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FIG. 4: (color online). Occupation probability of the dressed 
states \2a) (black upper traces) and \0a) (blue lower traces) 
(a). In (b) the probability for an atom to be in the untrapped 
dressed states is shown. The circles show exact numerical re- 
sults fSec. lllTAl . the dashed lines from the Magnus approxi- 
mation (Sec. IlllTjl . and the solid lines from the second RWA 
(Sec. If 11 CI . The calculation has been done at the rf trap 
bottom (x m i„,0, Zmin) for A = 1.05fii and 0,2 = 0.05f2i. 



analysis of Eqs. Ql-JHJ i s more complicated and conse- 
quently we evaluate these equations numerically. Some 
results for the time evolution of the populations of the 
adiabatic states \Qa) and \2a) are given in Fig. 0] We 
show both Magnus approximation (dashed lines) and the 
exact numerical solution (circles) for comparison. As can 
be seen, if we constrain the dressed spin dynamics to half 
of the first period, we obtain a very good agreement be- 
tween the first order Magnus series and the exact numer- 
ical results. After this time we see dephasing between 
the two evolutions and, even more dramatic, an impor- 
tant disagreement in the amplitude of the observed os- 
cillations. These two behaviours are somehow expected 
because at time instants very far from t — the con- 
tribution of the next order terms in the Magnus series 
becomes more relevant 26] . In Sec. IIII CI we can find a 
better approximation (using a second RWA) to the nu- 
merical solution which is shown in Fig. 0] with solid lines. 

Another test for the validity range of the first order 
Magnus series approximation is presented in Fig. [SI where 
the population of the dressed state \2a) and the proba- 
bility of leaving the rf trap are shown. In this figure 
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FIG. 5: (color online). Population of the dressed state \2a) 
versus position in the x (a) and z (c) directions at ilit = 20n. 
The probability for an atom to be in the untrapped dressed 
states is shown in (b) for the x and in (d) for the z spatial 
dependence (solid line - Magnus solution, circles and squares 
- numerical solution). As before, A = 1.05fh and O.^ = 
0.05fh. The inset shows a zoom of the resonances around the 
avoided level crossing at z m in ~ —0.19 mm. The solid line in 
the inset is just to guide the eyes. 
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the results of the numerical treatment and those from 
the Magnus series are represented by the points and the 
solid lines, respectively. To obtain the results in Fig. 03 
the initial adiabatic spin state has evolved during a time 
interval approximately equal to one half of the first oscil- 
lation period observed in Fig.0J in short, up to fiii = 20ir 
that is t = tt/0,2- Besides the good agreement that both 
methods show in the regions of less interest for us, we can 
notice the four resonances in the z dependence. The inner 
peaks are rather smaller than the outer ones, indicating 
that we have a position-dependent resonant coupling. 

C. Effective Hamiltonian from a second rotating 
wave approximation 

The first order Magnus approximation predicts well 
the location of the resonances and their spin rotation 

I 



where all the parameters appearing in it have already 
been introduced. We note that this equation is valid for 
any value of A, including those close to Q±, i.e., when w\ 
and u>2 are not so different from each other. Now we ap- 
ply a second rotating wave approximation, which is valid 
provided that the "detuning" A — Q and the maximum 
coupling SI2 are much less than the "carrier frequency" 
A. We finally obtain from Eq. I|1U|I the effective Hamil- 
tonian 

HA = -(A~n)Fo + ^[l + cos(9)]F ±e . (11) 

This last equation provides a new compact and power- 
ful description of the spin dynamics in the dressed trap 
in the presence of a second rf field. As an example we 
have shown in Fig. 0] the spin evolution (solid lines) pre- 
dicted with the effective Hamiltonian which is in excellent 
agreement with the exact numerical calculations. 

The form of the effective Hamiltonian (|llf> is com- 
pletely equivalent to that of Ha in Eq. I|A7[) . If we look, 
for instance, at the vectorial representation of the spin in 
the case of a single rf field (Fig.QJ, then in the presence 
of the second rf field one gets the picture shown in Fig. [HI 
where now Fsa — cos(9A)Fg+sin(9A)F±g plays a similar 
role as Fg did before. The angle 9a is then 

W^- ^^f , whh fl A e [0,,] . (12) 

We can view the resulting precession as a second dress- 
ing of the dressed spin states [l^. O ne can obtain the 



frequency. However, it fails to describe correctly the dy- 
namics away from the resonance points A = ±f2. To 
tackle this problem, we used a different approach and 
derived more generally applicable analytical expressions. 
The approach is based on a second rotating wave approx- 
imation, performed on Eq. JQ) and expressed through the 
rotation ~H(t) = KaH^U^ where TZ A = exp(iAtF e /h). 
This transformation leads to the time-dependent Hamil- 
tonian H (t) given by 



(10) 

I 

doubly-dressed states by diagonalizing the Hamiltonian 
(|ll|l . The corresponding eigenenergies of these states are 
given by multiples of HQ, a where clearly 




The spin oscillation frequency observed in Fig. 0] is pre- 
cisely Qa- On resonance, the period of these Rabi os- 
cillations induced by the coupling in 111|) is then T = 
47r/5l2[l + cos(0)l in agreement with the prediction of 
the Magnus approximation. 

Looking at Fig.[fJ]we realize that 9a — tt/2 corresponds 
to resonant coupling with a maximum transition prob- 
ability to a state orthogonal to the eigenstates of Fg. 
This happens when fl = A [see Eq. (|12fl ]. As in the 
Magnus case, this condition is realized by the two values 
9q = arcsin(£li/A) and ir — 9q with 9q G [0, 7r/2]. Re- 
calling that when 9 = tt/2 we are exactly at the avoided 
level crossings in Fig. |2Ib), 9q indicates the location of 
the outer resonances OR while it — 9$ takes care of the 
inner ones IR. 

IV. EVAPORATION 

A. General remarks 

In the following two subsections we will consider 
two schemes for implementing the evaporative cooling. 
Firstly, we will look at a pulsed scheme in which a frac- 
tion of the atoms are spin flipped out of the trapped state 



H(t) = - [A-fi-fia sm{9)cos{At)]Fg+n 2 [cos(9) cos 2 (At) + sin 2 (At)] F ±e + 

n 2 [l - cos(0)] sin(At) cos(At)F Y , 
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FIG. 7: (color online). Distance between neighbouring inner 
and outer resonances around an avoided level crossing. The 
circles are from the numerical calculation and the solid line is 
obtained using Eq. 1141 with b' = 203 Gauss/cm. The points 
were obtained at t = n/Q,2 by measuring the distance between 
the resonances as seen in the inset of Fig. |5f_d) . 



FIG. 6: (color online). In the presence of the second rf field 
and applying a second RWA, the spin (black arrow) precesses 
around an axis obtained when: Firstly, u/i tilts Fz by 9 getting 
Fe- Secondly, u)2 tilts Fe by 8a getting Fqa- 



by a sudden switch off of the second rf field. Then, sec- 
ondly, we will examine a continuous scheme in which hot 
atoms are removed from the rf trap by adiabatically fol- 
lowing a doubly-dressed state. In both these schemes the 
hot atoms that are going to be removed have to reach the 
resonances at 0$ or 7r — 6q. If we can neglect the grav- 
itational potential, the energy of these atoms should be 
larger than about FH(A — f2i) with respect to the bot- 
tom of the dressed rf trap. This approximation, valid 
for relatively small , can be refined using Eq. (|13|l and 
taking gravity into account. As an example, for our typi- 
cal experimental setup ^lj an d ^2 = 0.05Oi, this energy 
is equivalent to temperatures of 0.21, 6.1 and 11 /iK for 
A equal to 1.05Oi, 1.25f2i and 1.4f7i , respectively. 

These limiting energies imply, of course, that the atom 
cloud will have a finite size determined by the location of 
the inner and outer resonances. In Fig. 0we investigate 
the distance Azn ea between the neighbouring inner and 
outer resonances as a function of A. As expected, the 
distance between the resonances goes to zero when A 
is reduced. In fact, since ^2 A, and if we assume 
a constant magnetic field gradient we can derive the 
approximate form of Azu es from A = £1 as 



Az Res = 2 



V^ 2 - 



a 



(14) 



where a = gF^Bb' /h. 

The last point we would like to consider here concerns 
the inner resonance observed in the z direction. One pos- 
itive aspect about this resonance is that when A is such 



that both the inner and outer resonances are close to the 
bottom of the dressed rf trap, the atomic cloud trapped 
in the adiabatic state \2a) can be evaporated from both 
sides. However, the negative point is that some atoms are 
transferred by the inner resonance into the state | — 2 a) 
and trapped around z = 0. If these atoms come back to 
the region of the avoided level crossing, then they will be 
energetic enough that heating of the coldest atoms will 
take place via inter-atomic collisions. Note that for the rf 
dressed trap geometry discussed here gravity favours the 
evaporation through the outer resonance because a lower 
atomic energy is required than for the inner resonance. 



B. Pulsed evaporative cooling 

We first look at the pulsed evaporation scheme, which 
has as its main idea the extraction of hot atoms, in a con- 
trolled way, from the rf trap at the resonance locations. 
We can do this because at these locations we have large 
transition probabilities between the rf dressed states as 
seen in Figs. Efe), (f). These transitions have been al- 
ready analyzed, firstly, by using the Magnus approxima- 
tion (Sec. 1111 5j) and, secondly, by using the second RWA 
(Sec. IIII Cll . They have also been interpreted with the 
vector picture in Fig. H3 as rotations about F$a ■ Hence 
we can see that if at a given time instant t the precession 
axis of has zero projection onto Fg, then this state 

vector will be orthogonal to the initial dressed rf state 
and consequently, a transition has taken place. Clearly, 
the second rf field can transfer hot atoms out of the ini- 
tial trapped dressed state. The pulse has to be repeated 
several times during the trap oscillation period to ensure 
an efficient evaporation of all the atoms with sufficient 
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energy to reach the resonances. The pulsed evaporative 
cooling scheme requires that we can discriminate between 
hot and cold atoms by affecting as little as possible the 
atoms in the region of the rf trap centre. This implies 
a careful choice of the pulse duration and amplitude as 
will be discussed below. 

We already noted that in some situations it may 
not be desirable to evaporate via the inner resonance 
(Sec. llVl . One way we can avoid this resonance in 
the pulsed scheme is to carefully chose the time dura- 
tion of the pulse. For example, we can see in Fig. |3ff) 
that the depletion at the inner resonance (IR) takes place 
later compared to the outer resonance (OR). This hap- 
pens because, independently of the orientation of the dc 
magnetic field, the coupling with the Rabi frequency 
in lllfl is spatially inhomogeneous since 9 depends on r. 
To investigate this inhomogeneity for different values of 
A, we have chosen in Fig. |S|a particular value t = ir/fl2 
and explore the ratio of the probability of being in the 
untrapped states at the two resonances. This is shown 
in the figure as a function of the frequency difference A. 
The second RWA can predict this relative effectiveness 
of the inner and outer resonances. From Eq. 1)11(1 this 
ratio can be determined analytically for A = Q and it 
is plotted in Fig. [8] with a solid line. The calculation of 
the numerical results in Fig.|S]is done as follows. For a 
given value of A a figure similar to Fig. |SJd) is plotted. 
Then, the ratio of the inner peak height to the outer peak 
height is found. This is the coupling strength ratio we 
are interested in. Since there is an excellent agreement 
between the exact solution and the solid line in Fig. |H1 
we can state that indeed the second RWA works well in 
the parameter range we have explored. Notice that for 
A 3> Sli we recover the expected result that only the 
resonant coupling at the frequency ioi will occur between 
the bare states. For intermediate values of A, we can 
clearly see that the particular choice t = tt/^ for the 
pulse duration allows a good discrimination between the 
resonances IR and OR. In the final evaporation stage, 
where A > Oi, if we wish to limit the IR excitation it is 
necessary to adapt the pulse duration. 

As mentioned above we should avoid introducing tran- 
sitions at the dressed rf trap centre. Because of such 
transitions we can see in Figs.^a), (b), (e) and (f) that, 
even when the detuning for coupling adiabatic states is 
red or blue (A < f^i or A > f2i), the population of the 
initial trapped dressed state \2a) is in fact modulated at 
the centre of the rf trap (avoided level crossing). Such 
a modulation can produce unwanted heating or losses 
and, to study this process, we introduce the modulation 
depth. This quantity is defined as the contrast of the os- 
cillations presented in Fig. BJa) as the black points (nu- 
merical result). In Fig. we plot the modulation depth 
as a function of the Rabi frequency of the evaporation 
rf field, and for two different values of A. As expected, 
the modulation depth increases with since larger rf 
power is available for coupling the adiabatic states. Also 
not very surprising, this modulation is more pronounced 
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FIG. 8: (color online). Ratio between the coupling strengths 
at the locations of the inner and outer resonances. The exact 
numerical solution is represented by the points whereas the 
solid line is derived from the second RWA with A = fi. For 
A large compared to fii, the coupling at the inner resonances 
goes to zero. 



when A is such that ui 2 couples adiabatic states close to 
the position of the avoided level crossing. 




0.01 0.02 0.03 0.04 

Rabi frequency, Q. /Q. 



0.05 



FIG. 9: (color online). Modulation of the dressed state \2a) 
population at the trap centre for A = 1.05fii (black open 
circles) and 1.25f2i (blue open squares). The solid lines are 
calculated analytically from the second RWA treatment. 



Taking into account the result presented in Fig. we 
devised a strategy which affects as little as possible the 
coldest atoms, while doing the evaporation. The idea 
is to ramp and A simultaneously to preserve a fixed 
modulation depth at the rf trap centre [27j • In Fig. ^3 
such a ramp is presented, where we have allowed for a 3% 
modulation level of the coldest atoms population. Note 
that the reduction of O2 has to be taken into account for 
the optimization of the pulse duration. 

In fact, since Q2 and A can be controlled indepen- 
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potential V A (r) is given by 



1.0 1.1 1.2 1.3 1.4 

Frequency difference, A/Q 

FIG. 10: Optimized Rabi frequency as a function of the fre- 
quency difference A. In this optimization, the population 
modulation at the rf trap centre is limited to 3%. 



dently, we can make individual ramps for each one of 
these parameters and manipulate independently the po- 
sition and strength of the evaporation resonances. Tak- 
ing into account that there is a finite modulation of the 
atomic population at the rf trap centre, even for A very 
different from Cl±, we propose to use a limited number of 
rf pulses to cool down the sample. For example, a pulsed 
evaporative cooling scheme has been developed demon- 
strating the achievement of the collisional regime in a 
beam of neutral atoms (2^]. However, this pulsed scheme 
uses pulses longer than the trap oscillation period unlike 
the scheme proposed here. 



C. Continuous evaporative cooling 

Our second scheme for the evaporative cooling of atoms 
in the dressed rf trap uses the second rf field in a contin- 
uous rather than pulsed mode. This situation is closer to 
the normal case of the evaporation of atoms in a dc mag- 
netic trap by a single rf field. Here, we simply use the sec- 
ond rf field as a tool to control the dressed rf trap depth. 
This trap depth corresponds to the energy required to 
reach the resonances (with respect to the dressed rf trap 
bottom) as considered in Sec. lIV Al 

In the usual case of continuous rf evaporation it is use- 
ful to look at the system using dressed states. In our sit- 
uation, the equivalent relevant basis is given by doubly- 
dressed states. As remarked earlier, these are found by 
diagonalizing H A , Eq. (|1 1|) . or we can write 



h a = n A F eA , 



(15) 



with the frequency fl A as given in Eq. I|13l) . This fre- 
quency determines the doubly-dressed potential hfl A (r) 
which we would like the atoms to follow for the evap- 
oration to proceed. If we include gravity the resulting 



V A (r) = -Fhtt A {v) + Mgz 



(16) 



An example of this potential is given in Fig. 1111 which 
shows the minimum at z — z m i n (i.e. at 9 = tt/2 or 
6 = 0), where the cold atoms will eventually collect, as 
well as the resonance regions IR and OR which form the 
"lips" of the doubly-dressed trap over which the hotter 
atoms must pass. (During the usual evaporation process 
the "lips" are subsequently lowered by ramping the rf 
frequency and the same procedure can be carried out 
here with the second rf field.) 



G -0.2 

> 
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FIG. 11: (color online). Doubly-dressed potential Va shown 
as a function of distance z — z m in from the dressed rf trap cen- 
tre at z m i n for ^2 = 0.05f2i and A = 1.25f2i as in Figs. OJe) 
and (f). The potential is only shown in the region of the 
singly dressed rf resonance where S — 0, with the associated 
inner and outer doubly-dressed resonances indicated with IR 
and OR as in Fig.^f). Only the relevant, i.e. lowest, state is 
shown for F — 2 and gravity was taken into account in this 
figure. The vertical unit corresponds to about 20 fj,K for the 
value of f2i given in Sec. ITD 

For this picture to be valid, we must have an adiabatic 
following of the vector Fg A as the atoms move about the 
trap. A general condition for this can be expressed as 



al9 A 



dt 



< n A 



(17) 



where 8 A is given in Eq. (|12fl and the right-hand mem- 
ber of the inequality l|17fl is seen to be just the energy 
separation between the doubly-dressed levels. In prac- 
tice this condition is rather easily satisfied for a singly 
dressed rf trap (see, for example, Ref. ^H), which is 
relevant for the bottom of the doubly-dressed trap as il- 
lustrated in Fig. El However, in order to also satisfy 
the adiabatic following condition (jl7(l at the resonances, 
we will find a new constraint that should not be too 
small. The analysis of adiabaticity is conveniently car- 
ried out in terms of a Landau-Zener parameter A such 
that Eq. (jT7|l implies that A> 1. To proceed, we use 
the definition l|12|) of 9 A in order to compute its time 



10 



derivative assuming that only 9 is time-dependent. Cal- 
culating the time derivative and evaluating (|17f) at the 
outer resonance location 9o, we find that in terms of the 
Landau-Zener parameter A = ^a/|^a| the adiabatic fol- 
lowing condition reads 



A = 



ni[l + cos(fl )]' 
4|6 |cos(6>o) 



> 1 



(18) 



where Sq is the time derivative of the detuning 5(r) evalu- 
ated at #0: and is proportional to an atom velocity. From 
the expression for A we can see that as soon as f2 2 is 
reduced, only the slow atoms will have their spin adia- 
batically following the axis Fg&, and in fact, for motion 
linearized over the resonance, the multi-state Landau- 
Zener analysis [2(J shows that the probability for an atom 
to be lost from the adiabatic state in a single pass is 
1 — [1 — exp(— ttA)] 2F . For the parameters of Fig. ^2 
the energy of the atoms E (measured from the trap bot- 
tom in temperature units) would have to have a value 
of E/ks =20 /zK for the probability of a non-adiabatic 
crossing to reach about 10 -5 for the outer resonance OR. 

It is clear that if we switch from the outer resonance 
OR (at 9 — 6*o ) to the inner resonance IR (at 9 = it — 9q) 
the adiabaticity condition l|18|) will be different. This 
is connected to both the spatially dependent coupling 
0,2 [l + cos(0)] /2, and the effect of gravity, which as seen 
in Fig.llllmake the IR "lip" higher. The weaker coupling 
at the inner resonance also means that the dynamics is 
less adiabatic at this point. In fact for fairly "hot" atoms 
one can contrive that the resonance OR is rather adia- 
batic whilst the resonance IR is rather diabatic. Together 
with the effect of gravity, this would mean that atoms 
can be evaporated out of the OR resonance whilst adi- 
abatic coupling through the IR resonance is prevented. 
As explained in Sec. IIV Al this can be useful to partially 
prevent the return of evaporated atoms to the resonance 
regions with subsequent collisions and heating. However, 
we note that if we want to reduce the final temperature 
by steadily reducing A, the couplings at the two reso- 
nances become more equal (as in the pulsed scheme) and 
less discrimination between the two evaporation zones is 
possible. 



V. CONCLUSIONS 

We have seen that we can employ a doubly-dressed 
basis [lsf for the analysis of a dressed rf trap with two 
rf fields provided the second rf field is sufficiently weak. 
Using this arrangement of fields we can create a scheme 
for the evaporative cooling of atoms in a singly dressed 
rf trap in a continuous mode. In contrast with the tra- 
ditional continuous forced evaporation scheme the idea 
of evaporative cooling based on the application of rf 
pulses with properly chosen durations and frequencies 
is also developed for a dressed rf trap. The duration 
of such pulses is essentially determined by the power of 



the rf field used for the evaporation, although the opti- 
mal pulse length changes from one location to another 
in the adiabatic trapping potential. When the trapping 
and evaporative cooling radio-frequencies are compara- 
ble, we found the evaporation to happen also via addi- 
tional resonances. Even if these resonances can enhance 
the evaporation process in both the pulsed and contin- 
uous schemes, we have to be careful that atoms evap- 
orated through these resonances do not come back and 
heat the cold atomic cloud. In this respect, the advantage 
of the pulsed scheme is to provide an additional control 
of the transition probability at the inner resonance via 
the pulse duration. However, the effect of gravity is to 
favour the evaporation via the outer resonance OR which 
is then located at a position of lower energy. Although 
the importance of this effect may depend on the partic- 
ular trapping geometry (pancake, ring traps), the evap- 
oration schemes proposed in this paper are quite general 
and should allow an efficient cooling of atoms directly in 
the various rf traps which have been proposed or real- 
ized m m d m m m . 

The main application of the work in this paper is to the 
evaporative cooling of atoms in a dressed trap. However, 
there are also applications concerning noise and stability 
of dressed rf traps where the carrier frequency uj\ is not 
perfectly monochromatic. This can be due, for instance, 
to contamination by stray fields. Then it is clear that 
if the frequency components next to the carrier are in 
the range of, let's say, A = l.Olfii - 1.05fii, they will 
empty the rf trap if they have enough power to do so or, 
in the best case, they will raise the temperature of the 
atoms. This means that if we look at the second rf source 
as a noise term (a sideband in the frequency spectrum), 
the results presented here can be used to estimate the 
damage it causes. 
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APPENDIX A: DERIVATION OF H(r,t) 

We will start the derivation of Eq. JTJ by considering 
that the total magnetic field experienced by the atoms 
consists of three contributions or terms. One coming 
from the inhomogeneous dc magnetic field Bd c (r) of the 
QUIC trap, a second term oscillating at the frequency 
toi, Bi(r, t) = £?oi e i cos(wit), associated with the adia- 
batic trapping potential, and a third term of frequency 
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U)2, B 2 (r, t) — Bq2&2 cos(w2i), responsible for the evap- 
orative cooling in the rf trap. Here, _B i and Bq 2 are 
the amplitudes of the fields whereas ei and e 2 are unit 
polarization vectors. Using these definitions and denot- 
ing by F the atomic angular momentum operator, the 
total Hamiltonian of our physical system Ht{v, t) can be 
approximated by 

H T {r,t) = gp^BF- [B dc (r) +Bi(r,i) + B 2 (r,t)]/ft , 

( A1 ) 

where gp and [Ib are the Lande factor and the Bohr 
magneton, respectively. If we assume that at every point 
r the direction of the dc magnetic field defines the local 
Z quantization axis then, for X polarized rf fields, the 
Eq. ijAlfl takes the form 



H T (r, t) = uj (r)F z + Vi(r, t) + V 2 {v, t) 



(A2) 



In <|A2|> Wo(r) = .gF/-t_B-B £ ; c (r)/?i is the Larmor precession 
frequency. The interaction Hamiltonian Vj (r , t) is defined 
by the expression 

V j (r,t)=Sl j (r)Fx{e i ^ t + e- i ^ t ) , j = l,2 (A3) 

being fl/(r) = gF^BBoj / (2H) the Rabi frequency. 

Given Hf(v, t), the dynamics of an atomic spin state 
\<p(r,t)) is governed by the Schrodinger equation 

ifi^^=Hr(r > *)|0(r,*)), (A4) 
which in the frame rotating at the frequency u>\ becomes 



d\ip(r,t)) r 

ih = - S(r)F z + 

at 

n\Vi{t)Tii + n\v 2 {t)n 1 ]\^{v, t)) . 



(A5) 



In Eq. (SP we have introduced the detun- 
ing 8(r) = u>i — wo(r), the rotating frame opera- 
tor IZi = cxp(— iu>\tFz/K), and the rotated state 



\ip{r,t)) = H\\(j)(r,t)}. If we consider the bare state ba- 
sis {| - 2), | - 1), |0), | + 2), | + 1)} of a spin-2 system, 
the matrix form of the rotated interaction Hamiltonians 
Tl\Vi(t)TZi and Tl\V2(t)lZi are respectively given by 



fii(r) 



and 



n 2 (r) 



iAt 



1 






j'AI 



3 -iAt 
2 e 







V 



3 iAt 
2 e 



3 „-iAt 
2 e 






1 

0/ 







3 iAt 
2 e 





-iAt 






,iAt 

o J 



where A = dj 2 — ^i and we have made use of the rotating 
wave approximation (RWA) by discarding the terms that 
oscillate at 2ui, 2uj2, and ui\ + u>2- In general, we find 
the dynamics of \ip(r,t)) in l|A5() to be described by the 
Hamiltonian 

H(r, t) = H A (r) + il 2 [F x cos(Ai) + F Y sin(Ai)] , (A6) 
where the adiabatic Hamiltonian -Ha(t) is defined as 

H A (r) = -S(r)F z +Q 1 F x . (A7) 

The time-independent Hamiltonian l|A7Jl can be 
rewritten as H A (r) = Q(r)Fg if we define O(r) = 
y/6{r) 2 + tti, cos(6») = -<5(r)/fi(r), em.{B) = fii/fi(r), 
and Fg — cos(0)Fz + sm(6)Fx- 
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